Sonic crystal lenses that obey Lensmaker's formula 
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This paper presents a theoretical study of the phenomenon of acoustic imaging by sonic crystals, 
which are made of two-dimensional regular arrays of rigid cylinders placed in parallel in air. The 
scattering of acoustic waves is computed using the standard multiple scattering theory, and the band 
structures are computed by the plane-wave expansion method. It is shown that properly arranged 
arrays not only can behave as acoustic lenses, but also the focusing effect can be well described by 
Lensmaker's formula. Possible applications are also discussed. 

PACS numbers: 43.20.-(-g, 43.90.-fv 
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When propagating through periodic structures, waves 
reveal a particularly important property. That is, the 
propagation of waves will be modulated by the periodic 
structures. As a result, the dispersion of waves will no 
longer behave as in a free space, and the so called band 
structures appear. 

The phenomenon of band structures was first inves- 
tigated for electronic systems and was put on a solid 
foundation in the context of Bloch's theoremP|. Since 
the central physics behind electronic band structures lies 
in the wave nature of electrons, by analogy an imme- 
diate incjuiry would naturally be whether such a phe- 
nomenon can be demonstrated in classical systems in- 
volving, for instance, acoustic and electromagnetic (EM) 
waves. Although this problem was early addressed by 
Brillouin then by Yariv and YehQ, this seemingly 
straightforward question has not received serious atten- 
tion, until late 1980s. In 1987, Yablonovitch Q and John 
Pi proposed that such a band structure is indeed possi- 
ble when EM waves travel in periodic dielectric media. 
Since then, the study of EM waves in periodic struc- 
tures has been booming, eventually leading to the ex- 
ploration of acoustic waves in periodic structures (e. g. 
Refs. |S0,llllESIi3)- This has induced the establish- 
ment of the field of sonic crystals (SCs). 

The exciting phenomenon of band structures in sonic 
crystals allows for many possible applications. It has 
been recognized that SCs could be used as sound shields 
and acoustic filters 0, S El El El E3 ■ These appli- 
cations mostly rely on the existence of complete sonic 
bandgaps in which acoustic waves are prohibited from 
transmission in all directions. Thus most of early studies 
have been focused on the formation of bandgaps. 

Recently, interest in the low frequency region, in which 
the dispersion relation is more or less linear, has started. 
Since the wavelength in this region is very large compared 
to the lattice constant, the wave sees the media as if it 
were homogeneous. Consequently, it has been proposed 
that refractive devices could be developed to converge 
acoustic waves by SCs, that is the sonic crystal lenses. 
A necessary condition to be satisfied to construct a sonic 



crystal lens is that the acoustic impedance contrast be- 
tween the SC and the medium should not be large; oth- 
erwise acoustic waves will be mostly reflected. Once this 
condition is satisfied, the converging lens can be either 
convex or concave depending on whether the sound speed 
in the SCs is smaller or greater than that in the medium. 

Recently, it was experimentally demonstrated ^| that 
proper SCs can indeed make refractive devices such as 
Fabry-Perot interferometers and acoustically converging 
lenses. This interesting observation has stimulated fur- 
ther explorations of acoustic imaging by SCS|17, 18, 1.9|, 
and has been also extended to the optical imaging by 
photonic crystal lensesjUl^. In light of these novel de- 
velopments, in this paper we wish to explore further the 
properties of the sonic crystal made lenses. 

This paper presents a first attempt to quantify the 
imaging effect of sonic crystals. Here, we carry out nu- 
merical simulations on the focusing of acoustic waves by 
SCs. To comply with the current experiments [l^ so that 
experimental verifications could be readily done, we will 
consider the SCs made of arrays of rigid cylinders in air. 
It will become clear that the present approach can be 
readily extended to other SCs and also to photonic crys- 
tals. The results will demonstrate that SCs can make 
refractive lenses which nicely obey Lensmaker's formula. 
It is also shown that the refraction dominates the focus- 
ing properties. 

Before continuing, it is worth noting here that pursu- 
ing acoustic lenses can be dated back to 70s. For ex- 
ample, ultrasonic lenses made from bulk materials were 
proposed by Beaver et al. [Ill and Szilard et al. [23j |. 
However, these lenses are quite different from the sonic 
crystal lenses. The essential principles that govern the 
focusing phenomenon by sonic crystals are (1) the wave 
propagation are subdued by the crystal structure, lead- 
ing to the refraction of waves commonly seen at the in- 
terface of two different materials; yet (2) the frequencies 
of travelling waves are still within the passing band, and 
therefore no propagation inhibition will incur, allowing 
efficient transmission through the crvstalfl^. 

Consider N straight cylinders located at fl with i = 
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1, 2, ...A'' to form regular arrays. An acoustic line source 
transmitting monochromatic waves is placed at fl- The 
scattered wave from each cylinder is a linear response 
to the incident waves which are composed of the direct 
wave from the source and the multiply scattered waves 
from other cylinders, and subsequently contributes to the 
total waves, forming a self-consistent multiple scattering 
scheme. The final wave reaching a receiver located at 
TV is the sum of the direct wave from the source and 
the scattered waves from all the cylinders. Such a scat- 
tering process has been formulated exactly by Twersky 
[i^l , and was detailed in Refs 0, 0| . For regular arrays 
of cylinders, the band structures are computed by the 
plane-wave method to obtain the acoustic phase speed 
inside the arrays. 

In the following computation, the rigid cylinders of 
radius a are arranged to form triangular arrays in air 
with lattice constant d. The following parameters from 
the experiment JI^ are used in the simulations: (1) the 
radius of the cylindrical rods is 2.0cm; (2) the lat- 
tice constant is 6.35cm; (3) The filling factor, defined 
as the area occupation of the cylinders per unit area 
/s = (7r/2\/3)(d/a)2, is equal to 0.3598; (4) the sound 
speed in air is Va = 351m/s. Additionally, the acoustic 
transmission through the cylinder arrays is normalized 
as T = p/po, and the incident wave in this entire project 
is along the TX direction. We note here that the rigid 
cylinders can be any solid rods. As shown in Ref. any 
material whose acoustic impedance with respect to the 
air exceeds roughly 10 can be used as the composition 
material for the rods. 




FIG. 1: The band structures computed by the plane wave 
expansion method for triangular lattice of the rigid cylinders 
in air. 



Figure n presents the band structure result for the tri- 
angular array of rigid cylinders. We see that the disper- 
sion in the low frequency region is linear. Our subsequent 
calculations will be restricted to this region only. We cal- 
culate the sound phase speed from the linear dispersive 
curve as Vp,c = oj/K = 295m/s. Then the refraction 



index is calculated as n = Va/vp^c = 1.1898. We note 
that when waves propagate in structure media, it is usu- 
ally necessary to distinguish the group velocity, roughly 
characterizing the energy flow, and the phase velocity. In 
the linear dispersion regions, the two velocities are nearly 
the same. We have also computed the phase velocity in 
other directions, yielding the consistent values. These re- 
sults imply that the cylinder arrays may be regarded as 
an effective refractive medium, making it possible that 
SCs can be used as the refractive acoustic devices. This 
will be supported by the following simulation of the sonic 
crystal lenses. 

A conceptual layout of the acoustic lensing system 
made from a sonic crystal is shown by Fig. ^ The ge- 
ometrical factors, important to the following discussion, 
are clearly indicated in the figure. The source is place 
at a distance So from the lens. The image is at a dis- 
tance of Si. The cylinders, represented by the black dots, 
are placed insider the lenticular area. Inside the area, 
the cylinders are arranged to form the triangular lattice 
with the parameters given above. Fig. [factually shows a 
double two-dimensional convex lens of a refractive sonic 
crystal whose surfaces have radii of curvature Ri and i?2- 
We note that a three-dimensional lens can be fabricated 
by aligning two perpendicular two-dimensional lenses. 
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FIG. 2: The conceptual layout of an acoustic convex lens, in 
line with the optical lens depicted in Fig. 5.16 of Hecht|25| 



We first consider the acoustic transmission through the 
lenticular structure of SCs. The results are presented in 
Fig. 01 The scenario is illustrated by Fig. OJa). The two 
dimensional spatial distribution of the normalized trans- 
mitted intensity is shown in Fig. Efb) . Here the x axis is 
a horizontal line towards right, and the y axis is placed 
vertically upward. We see that the focusing of the trans- 
mitted wave is evident. Since the exact location of the 
maximum intensity point is not clear from Fig. |3fb), we 
have computed the variation of the field intensity along 
the X axis with y = 0, and also along the y axis with 
X ~ 20m. The results are shown in Figs. Etc) and|Sl[d). 
Fig. OJc) shows that there is a peak at a; ~ 20m. The 
intensity variation along the y axis shows a significant 
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peak at y = 0. Here we observe that the intensity is bet- 
ter confined along the y axis than along the x axis. Some 
relevant parameters are listed in the figure caption. 
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FIG. 3: (a) The line source denoted by S and the lenticular 
arrangement of rigid cylinders; (b) two dimensional spatial 
distribution of the transmitted intensity on the right 

hand side of the crystal; (c) the variation of the intensity 
along the x axis at y=0 and (d) the variation of the intensity 
along the y axis at x=20m. Parameters: (1) the source is 
place at 21m away from the center of the crystal; (2) the 
origin is placed at the rightmost point of the crystal; (3) the 
frequency is taken as 1500Hz, and the incident wave is along 
the TX direction; (4) the radii of the curvatures of the lens 
are chosen to be equal, i.e. R\ = R2 = 4m. The thickness 
of the lens is about 80 cm or 12 lattice constants, while the 
height amounts to 340 cm or 53 lattice constants; therefore it 
is a thin lens. The total number of scatterers is 546. 



where n is the refraction index. 

We have verified that the acoustic imaging by the SCs 
discussed above does obey Lensmaker's formula. First 
from the band structure, we obtain the phase speed and 
subsequently the refraction index as the ratio between the 
speeds in the air and in the SCs. Then the focal length 
can be calculated from Eq. (0) with known Ri and i?2. 
This will give the theoretical value for the focal length. 
The focal length can also be obtained by simulation. We 
vary the location of the source and then find the image 
position. Taking these two values into Eq. we can 
get a simulated focal length for each source location. If 
the device indeed behave as a lens, the simulated focal 
lengths at various source locations should be consistent. 
Furthermore, if the lens is refractive, the simulated values 
should also be consistent with the theoretical value. We 
found these claims are indeed valid for SCs. 
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FIG. 4: (a) The focal lengths versus the object lengths; (b) 
image distance versus the object length So- Both radii of the 
curvatures are taken as 4m. The theoretical results are shown 
by solid lines. 



Fig. El clearly demonstrate the acoustic converging ef- 
fects by SCs. Since the lenticular shape is thin, the fo- 
cusing effect may be quantified. According to the refrac- 
tion theory, it is possible to relate the object distance or 
length So, the image distance Si, and the focal length / 
by the thin formula , 
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If the medium can be regarded as a refractive medium, 
the focal length will be related to the radii of the curva- 
tures of the convex surfaces bv[25| 
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Eqs. and l(2Jl lead to Lensmaker's formula 
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FIG. 5: The focal lengths versus the radius of curvature at fre- 
quency 1500Hz. The crosses represent the simulation results, 
and the straight line refers to the theoretical estimates. 



As an example, we consider the lenticular shape of ar- 
rays of cylinders described in the figure caption of Fig. |31 
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The results are shown in Figure 0| Here, both results 
from the theoretical and simulations are presented. In 
(a), the focal lengths estimated from the simulation for 
various object distances are shown by the crosses. As the 
comparison, the theoretical focal length calculated from 
Eq. ||2Jl is also shown. In (b), the image distance versus 
the object distance is plotted. Here we see that (1) from 
(a) the simulated focal lengths are consistent with each 
other, and also consistent with the theoretical value; (2) 
from (b) the relation between the image distance and the 
object distance can be well described by Eq. Q and also 
agrees remarkably well with the theoretical curve. 

We have further verified the above agreement by vary- 
ing the radii of the curvatures of the convex surfaces. The 
comparison with the theoretical prediction is shown in 
Fig. El In the simulation, the thickness of the lens is kept 
unchanged at about 80 cm, when the radii are varied. We 
also take that both the radii of the curved surfaces are 
the same; so the theoretical focal length is / = 
The simulated value for the focal length is obtained as 
follows. For a fixed radius of the curvature, we vary the 
object distance and obtain the corresponding image dis- 
tance. Then the focal length can be obtained through 
Eq. 1^ for each object location. The focal length will be 
averaged over all the locations. Then we change to other 
values for the radius of the curvature. In this way we 
obtain the simulated focal lengths as a functions of the 
radius. The theoretical results are obtained from Eq. (0) 
with the refraction index obtained from the band struc- 
tures. Figure O shows that the estimated focal lengths 
are in excellent agreement with the theoretical predic- 
tion. Here we would also like to note that the agree- 
ment will deteriorate when the size of lenses is too small. 
In fact, when the size is too small, the diffraction effect 
will become dominant and therefore the focusing can no 
longer be described by the refraction process. 

In summary, here we demonstrate that SCs not only 
can make acoustic focusing devices, but also the focusing 
behavior can be well described by Lensmaker's formula. 
This finding will help design of actual acoustic refrac- 
tive devices, such as focusing audio speaker systems or 
sonars for large scale oceanographical probing purposes 
[26^1. Underwater sonars with sonic crystal lenses may 
improve the detection of fish 27] and other marine organ- 
isms 28|. Although the present work has been focused 
upon the acoustic cases and on the rigid cylinders, it is 
easy to see that the ideas can be readily applied to other 
SCs, and as well as the photonic crystals. 
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